Coulomb Drag between Quantum Wires 
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We study Coulomb drag in a pair of parallel one-dimensional electron systems within the frame- 
work of the Tomanaga-Luttinger model. We find that Coulomb coupling has a much stronger effect 
on one dimensional wires than on two-dimensional layers: At zero temperature the trans-resistivity 
diverges, due to the formation of locked charge density waves. At temperature well above a cross- 
over temperature T* the trans-resistivity follows a power law p oc T x , where the interaction-strength 
dependent exponent x is determined by the Luttinger Liquid parameter K c - of the relative charge 
mode. At temperature below T* relative charge displacements are enabled by solitonic excitations, 
reflected by an exponential temperature dependence. The cross-over temperature T* depends sen- 
sitively on the wire width, inter-wire distance, Fermi wavelength and the effective Bohr radius. 
For wire distances d S> kp 1 it is exponentially suppressed with T/Ef ~ exp[— dfei?/(l — K c -)\- 
The behavior changes drastically if each of the two wires develop spin gaps. In this case we find 
that the trans-resistivity vanishes at zero temperature. We discuss our results in view of possible 
experimental realizations in GaAs-AlGaAs semiconductor structures. 
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I. INTRODUCTION 

Measurements of Coulomb drag trans-resistivity be- 
tween two coupled low dimensional electronic systems are 
a powerful probe of scattering and correlations between 
electrons JlJ. In a measurement of the trans-resistivity 
Pb a current I\ is driven in one (the " active" ) of the 
systems, while no current is allowed to flow in the other 
system (the "passive" system). The Coulomb interaction 
between electrons in the two systems transfers momen- 
tum from the active system to the passive one, where a 
voltage drop V% develops. The ratio —V2/I1 is the trans- 
resistance, which is related to the trans-resistivity by a 
geometric factor. 

In weakly coupled two dimensional systems, at least at 
zero magnetic field, the trans-resistivity is usually pro- 
portional to the electron-electron momentum relaxation 
time, and is therefore proportional to T 2 , with T being 
the temperature. As explained by Fermi liquid theory, 
the T 2 behavior holds also in the presence of electron- 
electron interaction within each of the two coupled sys- 
tems. 

In one dimensional systems, which are presently re- 
alized by organic quasi-lD-metals, carbon nanotubes, 
edge states of quantum hall liquids, and ID semicon- 
ducting structures, electron-electron interaction is be- 
lieved to invalidate the Fermi liquid picture, and gen- 
erate a different state, described approximately by the 
Tomanaga-Luttinger (TL) theory (for reviews see 
e.g. Since electronic correlations in this state are 

stronger than in a Fermi liquid, it is interesting to exam- 



ine the Coulomb drag trans-resistivity between two such 
systems. 

In this paper we study theoretically Coulomb drag 
between two identical parallel one dimensional wires at 
close proximity For perfectly clean wires, as assumed 
here, the current flowing in the active wire generates volt- 
ages on the two wires, which, due to Galilean invariance, 
are equal in magnitude and opposite in sign. The trans- 
resistivity pd = — j^j- (with L the length of each wire) is 
then also the intrinsic resistivity (not including the con- 
tact resistance) of the active wire. Thus, we occasionally 
refer to pd as the "resistivity". Note, however, that this 
resistivity does not influence a symmetric flow of current 
in the two wires. 

Drag between ID electron systems was considered 
earlier by several authors. Hu and Flensberg M and 
more recently Raichev and Vasilopoulos M investigated 
the problem in the absence of electron correlations 
(apart from screening effects) within Fermi liquid the- 
ory. Tanatar ||To|| studied the same problem in the pres- 
ence of disorder, and Coulomb drag of Luttinger liquids 
with a point-like interaction region was considered by 
Flensberg [jllj, and Komnik and Egger flr2| . In a recent 
work by Nazarov and Averin |Ts) ] ID systems of spin- 
less electrons are treated as independent Luttinger liq- 
uids with coupling limited to inter-wire backscattering 
(Afc w 2kp , where fcp- is the Fermi wavevector in the two 
wires). The present work treats both intra- and inter- 
wire electron-electron interaction on equal footing. We 
find that although drag takes place primarily through 
2kF scattering, the small momentum component of the 
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inter-wire interaction and spin-density interactions affect 
it in a crucial way. The problem under consideration is 
also closely related to the problem of a coupled double 
(or N— )chain ]14|Jl5|]. In case of the spin- full problem 
results of Lee et al. JL4[ are useful. 

The paper is organized as follows: In section [0] we 
define the problem and present the main results and the 



physical picture. Section III deals with two wires of spin- 
less electrons (throughout the paper we use "wire" as 
synonym for "ID electron system"). After formulating 



this problem in subsection III B , we analyze it by means 
of a renormalization grou p in ill C| . We then discuss the 
high tempe rature regime ( III D ), the crossover tempera- 
ture ( III E| ) , and the low temperature regime ( III F ) . In 
this analysis we employ the method developed recently 
by Nazarov and Averin [ fl3[ , and earlier results on impu- 
rity pinned charge density waves (CDWs) [ ^8|JlS|l . In sec- 
tion IV wc address a double-wire system of spin unpolar- 
ized electrons. We write the Hamiltonian and renormal- 



ization group equation in subsection IV A, analyze the 
fixed points in IV B , deal with weak interactions in IV C , 



and extend the discussion beyond that limit in subsection 



IV D . We then discuss the high and low temperature lim- 
its of the drag in this case in subsection [V E . In section |v| 
we estimate experimental values of the relevant param- 
eters for semiconducting wires [p0|]2l|. section VI con- 



20[|21|. section 
eludes with a summary. Some technical details are put 
into the appendices. In particular, in appendix^ we ex- 
amine the relation between the Nazarov-Averin method 
and earlier weak coupling calculations of Coulomb drag 



II. REVIEW OF THE MAIN RESULTS 

We consider two identical wires of diameter d, sepa- 
rated by a distance d. We denote by kp the Fermi wave 
vector in each wire, by vp the Fermi velocity in each wire. 
The strength of the Coulomb interaction is characterized 
by r s = r/aB, with r the mean (intra- wire) electron dis- 
tance and as the effective Bohr radius. The length of 
the two wires is L. 

We first consider two wires of spin-less electrons. Ex- 
perimentally, this system may be realized by applying a 
magnetic field parallel to the wires, which would polarize 
the electrons' spins without affecting their orbital mo- 
tion. This system is closely related to a single Luttinger 
liquid with a spin-degree of freedom, when the two spin- 
projections are identified with electrons in the two wires. 
Therefore, results obtained previously on the effect of 
backscattering in such systems jl6],[l7| can be used. 

We find that, for infinitely long wires, Coulomb cou- 
pling always leads to a diverging resistivity pp> as tem- 
perature goes to zero. The physical picture behind this 
effect is that at sufficiently low temperature the electrons 



in both wires form two inter- locked CDWs. Then a rela- 
tive charge displacement can be created only by overcom- 
ing a potential barrier. At zero temperature this cannot 
be done by an infinitesimal electric field, and leads to a 
non-linear trans-resistance. At finite T, below a crossover 
temperature T* (discussed below), the trans-resistance 
satisfies, 

p{T) ~po lT exp(E a /T) 

with E s ~ T* defined below. 

For short wires a qualitative different behavior ap- 
pears. Here, from time to time the CDW in the active 
wire slips as a whole relative to the CDW in the other 
wire. These instantaneous slips are a result of either ther- 
mal fluctuations or tunneling events. The latter leads to 
a non-diverging resistance at zero temperature, which is 
exponential in L. 

At temperatures well above T* the previous picture of 
inter-locked CDWs is no longer valid. In this case it is 
more appropriate to think of independent electrons in the 
active wire, which suffer from backscattering at the 2kp 
component of the potential generated by density fluctu- 
ations in the passive wire. A perturbative calculation 
yields in this case a resistivity 



'" ri = '*> AL, ( Jr 



x = 4JC r - 



(1) 



exhibiting a characteristic power law dependence on T. 
The coefficient po is of order hkp/e 2 , A denotes the di- 
mensionless inter-wire backscattering potential, and E a 
is of order the Fermi energy. The parameter K c _ is the 
TL parameter of the relative charge density sector (c— ). 
It is determined by the difference of the small momentum 
intra- and inter-wire couplings, and not by the intra-wire 
small momentum coupling, as assumed in p3[ . With 
vanishing small-momentum interaction, K c - approaches 
unity, and p(T) takes the linear temperature dependence 
of the drag resistance of independent ID electrons [||. 
In the presence of Coulomb interactions K c _ may be ei- 
ther larger or smaller than 1, depending on the inter- wire 
distance. 

The crossover temperature T* is a complicated func- 
tion of four length scales: the wire separation d, wire 
width d, effective Bohr radius as, and the mean (intra- 
wire) electron distance r = ir/kp. The first, d, controls 
A, the strength of the 2kp component of the inter-wire 
Coulomb interaction. For widely separated wires d>r, 
this component is exponentially small, A cx exp— Ikpd 
and consequently, 



E exp - 



kpd 

1 -K r - 



As djr — > oo, then, T* — > and the trans-resistivity fol- 
lows Eq. ([l]) in all practically relevant temperatures and 
length scales. The general trends are shown in table Q. 
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The maximum values of T* that can be expected in re- 
alistic experimental set-ups are of order T* ~ 0.01 x Eq. 
In case of small wire separation d <C r the crossover tem- 
perature T* is exponentially suppressed according to 



T* ~ E exp 



7T 3 c(k F ) 
r s kpd 



where e(fcp) is of order one and only logarithmically de- 
pends on kp. 

For the spin- full case the results are similar, as long 
as the spin sectors are not unstable towards a formation 
of an energy gap in their spectrum (spin gap). The re- 
sistivity diverges at zero temperature, and scales with 
temperature with an exponent 

x = 2AV - 1 

in the high temperature regime. A comparison with the 
previous exponent reveals that fluctuations in the neutral 
spin sector moderate the effect of the charged modes. 

The behavior changes drastically if the single wires de- 
velop spin gaps. In this case we find that the trans- 
resistivity vanishes at zero temperature. 



III. SPIN-LESS DOUBLE WIRE 



A. Notation 



We use the following notations: 

(k),a rw {k) : creation and annihilation operator of a 
right (r = +) or left (r = — ) moving fermion of mo- 
mentum k. The second index refers to the active 
(w = +) and passive (w = — ) wire. 



(k) : fermionic operators in 



real space representation. 

p rw : density of right/left moving fermions in wire w. 

n w = P-, w + P+.w ■ charge density of wire w. 

n c ± = ri- + n + : absolute (or symmetric, c+) and rela- 
tive (or antisymmetric, c— ) charge density of the 
double wire system. 

In general, the indices c+ and c— refer to quantities of 
the absolute and relative charge mode, respectively. We 
also use this convention for bosonic field, introduced in 
the next subsection. 

The notation we use for the coupling constants follows 
that of Voit's review Q. We use gi to denote intra- wire 
couplings, and gi to denote inter- wire couplings. The 
subscript i — 1 denotes 2k F scattering, i = 2, 4 denote 
small momentum scattering. 



B. The Hamiltonian in fermionic and bosonic 
representation 



In this section we consider two one dimensional wires 
of spin polarized electrons with equal densities. If the 
wire index is viewed as a z-component of an " iso-spin" , 
the Hamiltonian of the problem is that of an iso-spin- ^ 
ID system (but with interactions that are not SU(2)- 
symmetric) , and results obtained previously on the effect 
of backscattering in such systems |l6|,[nj can be used. 
The kinetic energy part of the Hamiltonian states 0] , 



H a = v F ^ rka\ w (k)a rw (fc) = — ^ ^ p rw (g)p rw (- g) ■ (2) 

rwk qrw 



The small-momentum transfer or forward scattering part 
of the electron-electron interaction is given by j| 

Hf = — ^ [S w ,w>92 + 8 w - w >g2]p+ w [q)p- w >{-q) 

ww'q 

If ^2 [fiw,w'94+5 w ,- w >g4]Prw(<l)prw'(-q) 
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rww'q 



and backscattering processes are described by 



Hb = / dx ^\ w (x)^_ w ,(x)^ +w ,(x)'i])- w {x) x 

ww' 

[S w ,w'gi + 6w,-w'gi\- 



It is convenient to switch to a standard boson 
representation by introducing bosonic fields (j) w (x) = 
-iirL- 1 J2 q q^e'^-^/^p+^q) + P - W {q)) with their 
conjugates n„(i) - L" 1 £ g e -^- a M/ 2 ( P+w (q) - 
P-w{q)) @- Throughout the paper we interpret the 
length a as the inverse Fermi wavevector 2n/kp. Physi- 
cally, the field <fi w (x) denotes the displacement of elec- 
trons in wire w, normalized in such a way that den- 
sity fluctuations 5n w (x) and current I w (x) are given 
by d x (f) w (x) = -%8n w {x) and d t (f> w (x) = %I w (x). 
The relation to the fermions ip rw is established by 
the Luther-Peschel transformation formula, ip r w(x) — 
(2na)~ 1 ^ 2 exp[ir(kFX — </) w (x)) + 8 w (x)], where 9{x) — 
n J* o dx>Il w (x') §. 

The total Hamiltonian H = Ho + Hf + Hf, sepa- 
rates into two decoupled parts, one describing abso- 
lute (symmetric) current and density, and one describ- 
ing relative (antisymmetric) current and density. The 
decoupling is obtained by means of the transformation, 



2-V2^ + ± <f,_), n c± = 2-V2 (n+ ± n _) @. In 



bosonic representation the two parts are 
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^±Jdx (^K c+ n 2 Il 2 c+ + ^-(d^ c+ f 
^ jdx (k c ^ 2 U 2 c _ + j^id^c-y 




dx cos 



(-52 - 92 + gi), 

(-.92 + 92 +9i), 



2nv c - 

v c ± =v F + (g A ± .9 4 )/27r 
u c± =v c± (l-U 2 ± )V 2 . 



(3) 
(4) 
(5) 

(6) 

(7) 

(8) 

(9) 
(10) 



(The signs in the definition of the small momentum cou- 
plings U c ± are chosen according to the conventions in 
Jl7|].) The fields <fi c ± describe fluctuations in the absolute 
and relative density via d x <p c ± = —2~2ir8n c ±. Accord- 
ingly, the relation to currents are dt4> c ± — 2~^irl±. 

The current in a drag experiment is a superposition 
of a symmetric and an anti-symmetric current. A sym- 
metric current flows without resistivity, due to Galilean 
invariance. Thus, the resistance results from the anti- 
symmetric part only, and is determined by the relative 
charge sector H c _ only. Formally this is manifested in the 
invariance of H c+ to spatially homogeneous charge dis- 
placements 4> c ± (x) — > 4> c ± (x) + ip c ± , an invariance which 
is absent in the backscattering potential (oc cos \/80 c _ ) 
in H c -. Consequently, we confine ourselves to the sine- 
Gordon type Hamiltonian H c - . 

The Hamiltonian i? c _ has two parameters, K c - (which 
may be expressed in terms of U c -, see Eq. (^)) and 
g~\. Our results are all independent of the sign of g±, 
which we take below to be positive. For a single wire 
with electrons of two spin directions (namely, for SU (2) 
symmetric interaction, g<x = §2) the sign of U c — is deter- 
mined by the sign of g±, i.e., by whether the interaction 
is repulsive or attractive. Here this only holds for two 
wires which are very close to one another. For larger 
inter-wire distances J7 C _ can become negative also for 
repulsive interaction, in particular Coulomb interaction: 
Then, g~2 — 92 — — — \nd/d, while the parameter g\ is 
independent of d. Thus, for large inter-wire distance d, 
the parameter U c - < 0. 

Below we confine ourselves to repulsive interaction, 
and discuss the case of wires at close proximity and that 
of well separated wires. 



C. Renormalization group analysis 

In this subsection we analyze the backscattering term 
of the Sine-Gordon Hamiltonian by means of a renormal- 



ization group (RG) analysis, and show that if the bare 
interactions are weak, and the electron-electron interac- 
tion potential decays with distance, the drag resistivity 
diverges at zero temperature. 

Let us first recall the main elements of an RG treat- 
ment for H c - (see e.g. pi). For small backscattering 
couplings, A = gi/2Tru c - -C 1, the RG equations are of 
the Kosterlitz-Thouless type (here we denote K = K c -) 



dX 
dx 



= (2-2K)X, 



dK 
Ax 



= -2YK A 



(11) 



where the parameter x = \nl/a is the logarithm of the 
renormalized momentum cut-off I -1 . The RG procedure 
starts with the bare couplings Ao , K at an initial momen- 
tum cut-off Iq 1 of order of a -1 and ends with renormal- 
ized couplings at a final cut-off l^ 1 = max{L _1 , T /u c +\. 
One method (due to Jose et al. p2[) to derive the 
RG equations is to expand the scale invariant corre- 
lation function ^4>o-(xi,r 1 ) e -2i4 >e -(x2,T2)^ in powers f 

the coupling A and to integrate out the large momen- 
tum degrees of freedom. After re-exponentiating the 
result one can then read off the RG equations (for de- 
tails see [^2[ and also ^3|). Within each RG step only 
the backscattering interaction g~i is treated perturba- 
tively, whereas the small momentum interaction parame- 
ters g2 — 91, <?4, 52, 9i, which determine the parameter K, 
are treated exactly. A different method, which leads to 
the same RG equations ( |ll| ) but may provide additional 
insights, is by means of mapping the Sine-Gordon Hamil- 
tonian onto the two-dimensional Coulomb gas problem 

@- _ 

The integral curves A (if) of the RG-flow (|11|) shown 
in Fig. (nh obey the differential equation 



dA 
dK 



K - 1 
K 2 \ ' 



and are of the form 



A(/v-. = v2(l+ln(tf/ifo)--i- + f 



1/2 



There are two types of stable fixed points to Eqs. 
(pd|). Fixed points of the first type are characterized 
by A = 0,K > 1, i.e., by zero drag at T = 0. The 
basin of attraction of these fixed points is the area below 
the separatrix X S (K) = \[2~{^-^- + log K) . Systems 
with bare couplings Ao, Kq inside this area scale towards 
weaker backscattering coupling A when temperature de- 
creases. Below we show that no realistic set of interaction 
parameters falls under this category. 

When the bare couplings are outside this region (i.e. 
Kq < 1 or Ao > A s (-Ko) ) renormalization to lower 
temperature drives the system into the strong coupling 
stable fixed point, where A — > 00 and K — > 0. In 
this case backscattering becomes dominant at sufficiently 
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low temperatures and freezes the phase </> c _ to a mini- 
mum position. Translated to the double-wire system this 
means that the charges adjust their relative displacement 
V20c- = <fi+ ~ 4>- m such a way that the 2kp inter-wire 
potential is minimal, i.e. the system forms two inter- 
locked CDWs. 

Under these conditions the drag is very strong, as 
pointed out by Nazarov and Averin The system's 

resistivity to a flow of unequal currents in the two wires 
becomes infinite, in the limit of zero temperature and 
infinite length. We elaborate on this subject in section 



(|IIFj). 

It is instructive to express the condition for weak and 
strong coupling in terms of the bare interaction param- 
eters gi,§i- For reasonably weak interaction, K ~ 1.6, 
the separatrix X S (K) is well approximated by X S (K) 
K — 1 Ki U c -. Within this approximation the requirement 
for strong coupling, X > X S (K ), states 



A > U°_. 
which is equivalent to 

{92-91) > (52-5i)- 



(12) 



(13) 



Rather than the absolute strength of the inter-wire in- 
teraction couplings, it is their difference gi — g~\ in com- 
parison to 02 — 9i which determines the zero temperature 
fixed point of the system. Taking and <?.; as the Fourier 
components of intra- and inter-wire interaction V(x) and 
V(x) at qi — 2kp and (72 = l/£ — * 0, the condition 
becomes 



dx (1 - cosgix) [V(x) ~ V(x)] > 0. (14) 



The first factor is non-negative. Then, a sufficient condi- 
tion for the l.h.s. to be positive is obviously V(x) > V(x) 
for all x, which is fulfilled by all monotonously decay- 
ing repulsive interactions potentials, in particular the 
Coulomb potential. Therefore, a Coulomb coupled dou- 
ble wire system of spin-less electrons should scale towards 
strong coupling, which implies a diverging zero tempera- 
ture drag in infinitely long wires. 

There are two types of initial values Ao, K$ which flow 
to the strong coupling fixed point. The first is defined 
by Kq < 1. As seen in Eqs. (|Tl|), for this case A varies 
monotonously as the temperature decreases. The sec- 
ond type is defined by K > 1 and Ao > X s (K n ). For 
this type, A does not varies monotonously. For relatively 
high temperature, A decreases. At the temperature at 
which K = 1 A starts increasing towards the strong cou- 
pling fixed point. As explained in the previous subsec- 
tion, well-separated wires {d — » 00) fall under the first 
category, while wires at very close proximity fall under 
the second. Assuming the bare interaction parameters 
gi,g~i to be small, for both types of initial conditions we 



may separate between a weak coupling, high tempera- 
ture, regime, where perturbation theory calculations can 
be carried out, and a low temperature, strong coupling 
regime. In the next subsections we calculate the drag re- 
sistivity for both regimes, and identify the temperature 
scale that separates the two. 



D. The high temperature regime 

We begin with the weakly coupled regime, in which we 
employ a method devised by Nazarov and Averin [ jl3| . 
In the limit of linear response (/ — > 0) this approach of 
calculating the drag resistivity is similar to the memory- 
function formalism of Zheng and MacDonald |25| , as it is 
shown in appendix [b]. The main difference is that in the 
present calculation only the backscattering component of 
the inter-wire interaction is treated perturbatively, while 
the small-momentum part is treated exactly. 

We consider a four-probe measurement with voltage 
probes at positions xq and xq + a on both wires (let 
w c _/T <C a <C L), and calculate the voltage drop 
eV w = (fA w (xo) — /j, w (xq + a)} j along wire w when a cur- 
rent / is driven only through the active wire (+). Using 
the relation 5eV w = n~ 1 5n w , where k~ 1 = d[i w jdn w is 
the inverse compressibility, we obtain, 

eV w = ~(8n w (x + a) - 6n w (x))j = — / (d x n w ) I . 

K K Jx 

Due to translational invariance, 

PV 1 7T 

= -(d x n w )j = ~ (dl^ + +wdl4> c -) r (15) 
a K V2k 

The thermo-dynamical averaging (. . .) j has to be re- 
stricted to states satisfying = I and (/_) = 0. 
Equivalently, but technically more convenient, one can 
use an ensemble of current-less states, and then perform 
a Galilean-transformation of the active wire such that a 
net current I results. In terms of the displacement fields 
<p s this means that </> + acquires a component growing 
linearly in time, 

4> + (x, t) — > 4> + (x, t) + Qt, 



or, translated into absolute and relative fields, 

fit 



r >c± — > <Pc± 



(16) 



The frequency VI is related to the current by Q = nl/e 

As expected, the transformation (llq) does not alter the 
absolute sector H c+ , and no symmetric voltage is induced 
by the current. It does affect the relative Hamiltonian 
H c - via the backscattering interaction, which becomes 



5 



dx 

Hi nt = ^Eq I COS 

ira 



and gives rise to a finite drag-voltage 



2Qt), 



a 



71 



(17) 



(18) 



For the following calculation it is advantageous |l3|] 
to make use of the equation of motion, df4> c - = 
— [H c _, [H c _, C _]], from which follows that under sta- 
tionary conditions ((<9 4 2 c _) = 0) 



K /SHj, 
u \ 5<f) 



(19) 



and then to perform a perturbative expansion of the 
r.h.s in the backscattering coupling A with respect to 

^0 = H c - — Hint- 

By standard methods we obtain in lowest non- 
vanishing order the resistivity 



d_Vu, 
dl I 



wp X 



1=0 



4K-3 



(20) 



where p ~ h/e a. Higher order terms p^ n > oc A scale 
with temperature as 



p 



[AK -4)n + l. 



Terms with odd powers of g\ vanish. 

The temperature scaling law p oc T 4K ~ 3 can be also 
derived from the RG treatment in the following way: In 
the absence of intra-wire interaction a simple calculation 
yields 

P (T) ~ p \ 2 T/E . 

The effect of intra-wire scattering can be accounted for 
by using a renormalized backscattering coupling constant 
A = A(T). One finds from the RG equations (O) that 
in the weak coupling limit A(T) = A (T/£ ) 2 (see 
appendix |A| Eq. (Al)), which inserted into the previous 
equation indeed gives p(T) oc Aq (T/E ) iK ~ 3 . 



E. The crossover temperature 

In this section we estimate the cross-over temperature 
below which the perturbative calculation of the previ- 
ous section ceases to hold, and inter-locking of the two 
charge density waves in the two wires becomes relevant. 
The cross-over temperature is that in which the coupling 
constant A becomes of order 1. 



For initial conditions Ao 1 and Kq < 1, correspond- 
ing to well separated wires, the solution to the RG equa- 
tions can be approximated by A(T) 
K w Kq. The cross-over is then 



Ao (fL ) 2 ~ 2Ka and 



p rl/(2-2Ka) 



(21) 



In appendix |a| it is shown that for ii > d, kp and 
for the Coulomb interaction, the coupling constant Ao ~ 
cxp—kpd, and thus T* is exponentially suppressed. 
More specifically, we find in this case 



T* r~j Eq exp 



kpd 
1-Kq 



(22) 



Note that with vanishing interaction strength, 1 — K c ^ s» 
— U c _ ps v s {\a.kpd + 7)7r -2 (see section^), and hence 
T* oc exp(— - d ), where the parameter r s = r/as — 
ire 2 /vpe describes the interaction strength. 

For initial conditions of the type Kq ~ 1 and Ao > 
X S (K), corresponding to wires at very close proximity, 
the RG equations can again be solved approximately (see 
Appendix (|A|) for details), and we find, 



E Q exp - 



7T 3 C(kp) 

kpd 



(23) 



where c(kp) is of order one and only logarithmically de- 
pends on kp. Thus, when 1/kpd 1, the temperature 
T* is again exponentially small. At first sight, this might 
be surprising. However, it becomes understandable by 
observing that for d <C r intra and inter wire coupling 
are almost equal. Hence the problem becomes close to 
that of a single wire with electrons of two spin states, 
(where gi — gi), which is known to renormalize to the 
marginal weak coupling fixed point. 



F. Low temperature limit 

In the strong coupling limit, with A ~ 1, and K <C 1, 
the Sine-Gordon Hamiltonian (^) describes two inter- 
locked charge density waves, and carries strong resem- 
blance to the problem of a pinned charge density wave, 
which has been studied quite extensively in the 1970s (we 
refer in particular to the work of Rice et al. |l8| and Maki 
0)- 

Classically, for large A the Hamiltonian H c ^ has an in- 
finite number of ground states with field configurations 
<j>(x) = 4>n — n(2N + l)/\/8- They are separated by a 
large energy barrier, such that at sufficiently low temper- 
atures the field fluctuates around one of the ground states 
<f>N- These fluctuations do not carry an anti-symmetric 
current. Rather, it is a transition from a ground state 
<f>jq to its neighbor 4>n±i which corresponds to an elec- 
tron transfer through the active wire from the left end 



G 



to the right end or vice versa. These transitions are car- 
ried out either by a soliton moving along the wire, or 
by a soliton-antisoliton pair which is created in the wire 
and is dissociated by an electric field. The width W and 
energy E s of a classical soliton are fl26J] 



E„ = 



8A 



-En 



We expect that the transport properties of short wires 
with L<lf are different from that of long wires, where 
L ^S> W, and therefore consider both regimes separately. 

In long wires electronic current is carried by means of 
soliton-antisoliton pairs. This mechanism was studied by 
Rice et al. [jl8| and Maki jl9) within a semi-classical ap- 
proximation. The semi-classical regime is that in which 
the classical soliton energy E s is much larger than the 
zero point energy of the fluctuations around it. In terms 
of our parameters this regime is defined by K -C 1. As 
found by Rice et al. p8[ , in this regime thermal creation 
of soliton-antisoliton pairs leads to a resistivity 



' E a T 
2XKE 2 ' 



(24) 



where I, the mean free path of (anti)solitons, is a phe- 
nomenological parameter. 

When K <C 1, the soliton energy is renormalized ac- 
cording to pTl] 



e f -'tr 



(25) 



Eqs. (24) and ( |25| ) define a crossover temperature T* = 
E r s separating the low and high temperature regimes. As 
one may expect, this temperature is approximately equal 
to the one obtained in Subsec. 
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The low temperature regime is characterized by a 
highly non linear drag, taking place at exceedingly low 
temperatures, where soliton-antisoliton pairs are gener- 
ated by quantum tunneling, rather than thermal activa- 
tion. In this regime, as found by Maki, 



e 2 16tt 2 u c 



e /e 



(26) 



where eo = E 2 /2eu c ^,e — V/L, and V is the voltage 
difference between the two ends of the wires. The cross- 
over temperature, at which thermally activated behavior 
changes to a tunneling dominated one, can be determined 
by equating the exponentials, which yields 



For short wires, where the soliton width W is large 
compared to the wire length L, the fields II, <j> may be 



approximated as constant in space, leading to an Hamil- 
tonian 



H = 



u c -L 
2vr 



n 2 K c ^U 2 + ^4 cos 
a 2 



E h 



(27) 



Here we introduced a frequency uj 2 = 8XKEq, a barrier 
energy Eb = -X—Eq, and canonical conjugated variables 
tp = \/8(f), p = LU/i/8. While this Hamiltonian resem- 
bles that of a pendulum, it differs from it in one impor- 
tant aspect: in the present case the states ip and (p ± 2?r 
are distinct physical states, since a shift of the phase by 
2ir corresponds to a transfer of one electron between the 
reservoirs the active wire is coupled to. Since different 
minima of the potential energy correspond to distinct 
states of the reservoirs, it is unlikely that the field <f> can 
be in a superposition of several such minima. We there- 
fore assume that the field is in one of the minima, and 
transport takes place by tunneling or thermal hopping 
between adjacent minima. 

A chemical potential difference <5/i between the reser- 
voirs at the ends of the active wire adds a linear potential 
to the Hamiltonian ( p7| ) 

To determine the current induced by Sfi we calculate the 
rates Ty r by which the phase tp hops between minima by 
thermal activation or tunneling. Then 

i = e (r r -r,). 

Thermal activated hopping dominates at temperatures 
T ^> uj with a rate IE8I 



r 



l/r 



2^ 



At low temperatures T<Cw the main contribution comes 
again from tunneling processes with rate 

T' l/r ~cje- s '/^. 

For S/j, <C Eb the action can be expressed as 



Sl/ r (Sn) = 



E b 
V2l, 

UJ UJ 



dtp\ 1 — cos ip ± — — 
■nE b 



where a ~ (Ka 2 /XL 2 ) 1 / 4 <§; 1 takes care of the finite 
width of the well ground state. The coefficient C is of 
order unity and only logarithmically dependent on a. 

Putting this together we end up with current-voltage 
relations 
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I « e -e- Eb/T swh(5p,/2T), for T > w 



J~ ewe- 4Bi>/a; sinh(C5At/w), for T < w. 
Accordingly, the linear resistance is 

R^^-e Eb/T , for T>u, 



for T < w. 



(28) 
(29) 

(30) 
(31) 



These expressions are valid for w ires much shorter than 
the soliton length W ~ a/V KX and voltages eV <C 
Notice that because of Ef, cx L/a the resistance increases 
exponentially with length L. 



IV. SPIN-FULL DOUBLE WIRE 

A. The Hamiltonian and renormalization group 
equations 

We now extend our analysis to include the spin degree 
of freedom of electrons in the two wires. Denoting the 
spin degree of freedom by an index s = ± in addition to 
r,w as above, and taking into account the SU(2) sym- 
metry of the Coulomb interaction, we have 

TTVp 



Hq — —jr- Prsw{<l)Prsw{ — q), 



qrws 



Hf = j ^ [dw,w>g2+fiw,-w'92]p+sw(c[)P-s'w>(-<l) 



qsws'w' 



+ — [S w ,w f 9A + &w~w'g4\prsw{q)Prs'w'(-q), 

rqsws' w' 

H b = ^ J dx [8 w>wl gx + 6 Wi - w >gi] x 

sws'w' 

^+sw i^-s'w' (x)4>+s'W {x)^-sw (x) ■ 

Following the same procedure as in the spin-less case, 
we first change to bosonic field variables 

i n e -iqx-a\q\/2 
<l>sw{x) = ~— > (P+sw(q) + p-sw(q)), 

L q 

n«,(a:) = e- iqx - alql/2 (P+su,(q) - P- sw (q)), 

q^O 

and then transform to new bosonic fields cj) c / s ± , IL C i s ± 
corresponding to absolute (+) and relative (-) charge (c) 
and spin (s) density. The transformation is defined by 

4>sw = 2^ c + + S ^ s + + W 0c- + SW(j> s -), 

and in the same way for n c / s ± . For the transformation of 
the backscattering Hamiltonian Hb we use the formula 



Vv 



'2na 



exp [ir(kpx — <j) sw ) + i6 



(32) 



where sw {x o ) = ir 1 U sw dx and the r/ rsw = 7y| stu are 
so-called Majorana fermions. They obey 

and ensure that the operators given by (|3^) follow 
fermionic commutation rules [Q. We obtain the follow- 
ing Hamiltonian (the tilde is omitted, j takes the values 
j = c, s, for charge and spin. I takes the values I = ±, 
for symmetric and anti-symmetric) 



H — 2^ Hji + Hb, 



(33) 

(34) 
(35) 



K jl = {l + U i i)V\l-U jl )- i '\ 

U c+ = ^—(-2g 2 - 2g 2 + 9l ), (36) 



2ttv, 
1 



c+ 



U c - = - (-2.g 2 + 2g 2 + 9l ), (37) 

27TV C _ 



U s ± = 



1 



-9i, 



2irv s ± 

v c ± = V F + (.94 ± 94) 

V s ± = V F , 

Ujl = Vjl (l-UW 2 , 



H b = + 



91 



f 1 „ / dx cos20 s+ cos2</> s 
dx cos20 c _ cos20 s 



91 



dx cos20 c _ cos20 s+ . 



(38) 

(39) 
(40) 
(41) 

(42) 
(43) 
(44) 



Again, the system has no resistance to a symmetric flow 
of current, due to Galilean invariance, and thus the 
charge symmetric part of the Hamiltonian H c+ decou- 
ples from the other parts. The remaining parts are all 
coupled by the backscattering Hamiltonian Hb, so that 
the spin modes, despite their charge neutrality, do affect 
the role of the Coulomb interaction. 

As in the spin polarized case, when the backscattering 
couplings g\ , g~\ scale to the strong coupling fixed point 
at T — 0, spin fields (j> s ± and the relative charge field 
4> c - freeze to minimum positions, which causes a diver- 
gent drag resistance at zero temperature. On the other 
hand, when A renormalizes to zero (while A can still scale 
to stronger couplings) , the drag resistance decreases with 
temperature and vanishes at T = 0. (The case of a de- 
creasing A and an increasing A does not exist.) 

The second order RG-equations for a double wire as 
described by the Hamiltonian above are 
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^ = (2-K s -K c - 2A)A, 

^ = (2-2if s )A-2A 2 , 
dx 



dK s 

dx 
dK c 

dx 



= -2(A 2 + A 2 )K 2 , 



= -4A 2 if 2 



(45) 
(46) 
(47) 
(48) 



simple physical interpretation, separating between 
the region where bare intra-wire backscattering is 
stronger than the bare inter- wire one and the region 
where the opposite is true. 



C. The small U limit 



where A = — <7i/27rat c _, A = gi/2nu s ±, K c = K c -, 
and K s = K s ±. These equations can be derived in 
the same manner as those in the spin polarized case 
by assuming scale invariancc of the correlation functions 



<« 



2i<jiji{x 1 ,T 1 ) -2i(f>ji(x 2 ,T2)^ 



) US T - ey are valid for ar ~ 



bitrary K s / C , but restricted to small |A|, |A| <C 1. (Similar 
RG equations appear in the context of single wall carbon 
nanotubes, where also two spin-degenerate channels are 
present |29|.) 



B. Fixed points for the renormalization group 
equations 



In case of small U c i s Eqs. fiq ) to ( [48[ ) reduce to the 
following set of equations, derived first by Gorkov and 
Dzyaloshinskii for the problem of coupled chains pfl,H , 



dx 



+ U C + 2A)A, 



^ = -2t/ s A~2A 2 , 
dx 



d£c 
dx 



= -4A 2 



(49) 



The following types of fixed points are found for the 
RG equations ©-(liF 



Fixed points where K s ^ 0. As evident from (47) 



for these fixed points A = A = 0, and thus we re- 
fer to them as the "weak coupling fixed points". 
They describe wires with vanishing backscattering 
interaction at low temperatures, and therefore van- 
ishing drag. As in the spin polarized case, for large 
values of K c , K s these fixed points are stable, while 
for small values they are unstable. 

• Fixed points for which K s — and K c ^ 0. As 
evident from (f48|), for these points A = 0, and 
then, from ([l6]), a stable fixed point exists only for 
|A| = oo. The former, A = 0, indicates vanishing 
drag, the latter, |A| = oo, means that the spin- 
modes are massive (spin gap). These fixed points 
have no analog in the spin-less case. We refer to 
them as the " spin-gap fixed points" . 

• Fixed points where K s = K c = 0. There are 
two stable fixed points in this plane. The first 
is A = — oo and A = +oo, and it obviously de- 
scribes two strongly coupled wires, with diverging 
zero temperature drag. We refer to it as the "strong 
coupling fixed point" . The second is a spin-gap 
fixed point as described above. There are also two 
unstable fixed points on that plane. The first is 
A = A = 0. It is repulsive in all directions of the 
(A, A) plane. The second is A = A = 1. It is attrac- 
tive in the direction (1,1). The border between the 
basin of attraction of the two stable fixed points is 
the diagonal A = A (see Fig. ||). This border has a 



(U c / S = J7 c -/s±)- Furthermore, in that limit A = U°, 
and d(X - Us)/dx = 0. Therefore X(x) = U s {x). With 
this, Eqs. (E^) reduce to three independent equations: 



| = -(^ + 3A)A, 



^ = -2(A 2 + A 2 ), 
dx 



d£c 
dx 



= -AX 1 



(50) 



These RG flow equations were analyzed in detail by Lee 
et al. G3] , who found that for initial values satisfying 



U C X < A 2 , 



(51) 



the system scales to a strong coupling fixed point A — > 
— oo, A — » oo. Rewritten in terms of gi...,gi... this condi- 
tion states 



fj2 



s Mi 
92 < x( — 

2 gi 



9i) + 0(glJv 2 F ) 



For gi > the r.h.s. is always larger than g~\ — <7i, such 
that a sufficient condition for this inequality to be fulfilled 
is g~2 — 92 < 9\~9\- This, however, is exactly the inequal- 
ity (|l3|) we had in the spin-less case. We find then again 
that for weak monotonously decaying repulsive electron- 
electron interaction (small positive A), a double wire with 
spin unpolarized electrons scales towards strong coupling, 
manifested by diverging trans-resistance at zero tempera- 
ture. The condition Xq = U®, valid in the small U regime, 
does not allow a flow to the spin-gap fixed point. 
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D. Beyond the small U limit — linear stability of the 
weak coupling fixed points 



As stated above, a small U analysis (i.e., weak 
electron-electron interaction) of the RG equations (|4q)- 
(fl8|) leads to the conclusion that A and therefore also 
the drag resistivity diverges at zero temperature, ac- 
companied by a large negative coupling A. It is un- 
likely that stronger electron-electron interaction would 
push the system towards the weak coupling fixed points, 
where both backscattering interactions become negligible 
at zero temperature. It is, however, conceivable that for 
some range of parameters strong interaction would lead 
the system towards the spin- gap fixed point, where A 
goes to large positive coupling, but A (and po) vanishes. 
For example, if in the absence of inter- wire coupling each 
of the wires is in a spin-gapped state with A ~ 1 (as 
happens for attractive intra-wire electron-electron inter- 
action) one may expect weak inter-wire coupling to leave 
the two wires in that state. Consequently, the wires 
would decouple (A — > 0) and the drag would vanish. 

In principle one may analyze this behavior by use of the 
RG equations (p|)-(p|) beyond the small U limit. The 
validity of such an analysis to the problem at hand is how- 
ever unclear, since at the level of initial values Ao > Ufj , 
and the RG equations are derived under the assumption 
of small A. Instead of doing this, we treat K c / S , A and A 
as independent parameters and confine ourselves to the 
case of small couplings A, A -C 1 (but arbitrary K c / S ), 
where the RG equations are valid. 

Let us consider the linear stability of the weak coupling 
fixed points (where A = A = 0) to turning on a small A, A. 
We find that the A = A = plane is split into four re- 
gions (see Fig. |^). In the first region, area 7, defined by 
K s + K c > 2 and K s > 1, the weak coupling fixed point 
is stable. In the region where K s > 1 but K s + K c < 2 
(which we call area 77, see Fig. ||), the weak coupling 
fixed points are linearly stable with respect to infinitesi- 
mal values of A, but unstable with respect to such values 
of A. We now argue that initial values in that area flow 
to the strong coupling fixed point (A — > -co, A — > +oo). 
Let us consider |Ao| <C Ao <C K s — 1 in this region. By 
Eq. (phf), the coupling constant A will scale to small val- 
ues, such that the instability of the A-mode (Eq. |45| ) 
becomes relevant . Then A keeps increasing until finally 
the — 2A 2 -term in Eq. (|t|) overcomes the A-stability and 
forces A to scale to negative values. 

Points lying below the K s — 1 line in area III or 777' 
are linearly unstable in A. Consequently, A increases and 
thereby decreases the coefficient in the equation for A. In 
the limit of Ao <S Ao this leads to a stabilization of the 
A-mode, even where it was initially unstable (area III'). 
Hence, points in areas III and III' scale to the spin-gap 
fixed point (A — > +oo, A = 0). 

To summarize, in the limit of |A| <§; A <C \K S — 1| the 



basins of attraction of the weak and strong coupling fixed 
points correspond to areas I and II in Fig |j| respectively. 
Further, points of area III /III', which characterize sys- 
tems with tendency to (single-wire) spin gap phase, in- 
deed scale towards the (double-wire) spin-gap fixed point 
with vanishing A (zero drag). 

The linear stability analysis we carried out is expected 
to give the right structure of Fig. (||) , but not the precise 
borders between the areas. Near the borders terms linear 
in A and A may have very small prefactors, which make 
the role of the quadratic terms significant. An example 
to that is the region K s « 1 + A, studied by Lee et at, 
where in the right hand side of ([h3|) both terms become 
comparable. 

When Ao is comparable to \K S — 1| (but still Ao -C 1) 
the condition for A scaling towards weaker couplings is 
K s — 1 > Ao (see the discussion in Sect. III). Then, 
the border between areas III /III' and II /I is raised to 
K s = 1 + Ao- The effect of larger Ao is therefore a re- 
duction of area 77, i.e. also in this sense a larger Ao is in 
favor of wire decoupling. 



E. High and low temperature limits of the drag 

A perturbative calculation of the drag-voltage in the 
weak coupling regime, applicable for widely separated 
wires in the high temperature regime, can be developed 
along the same lines as before. We have 

eV w = — {dl(j) c+ +wdlcj) c -) r 

The transformation to current carrying states is again 
given by (|l6|), such that 



eV w = w — (dl4> c J) I 



4dK - 



7TCX 



— w — -A(sin(20 c _ + V2flt) (cos 20 s ^ + cos2(j) s+ ) 



Hereby, we used the Eq. of motion for </> c _ and the sta- 
tionarity condition (<9 2 c _) = 0. A perturbative expan- 
sion in the couplings A, A and U s ~ A then leads to 



Pw = wpo A 2 



2K.-1 



(52) 



with p ~ h/e 2 a. Leading higher order terms are of order 
A 2 x 0(\,U S ). 

For wire distances d ^s> kp 1 the smallness of A cx 
e -2k F d j g g uaran teed. Further, in this case the spin- 
sector couplings A, U s flow to the single wire Luttinger 
Liquid fixed point A* = U* = 0, as long as the temper- 
ature is still large compared to 7*. Hence, also effective 
couplings A w U s are small, wherefore the perturbative 
result (p2h is applicable. 
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The cross-over temperature T* , separating between the 
high and low temperature regimes, is exponentially sup- 
pressed for J 3> kp l with 



T*~£ exp - 



2k F d 
l-K r 



(53) 



as shown in appendix |A| 

When the system flows towards the strong coupling 
fixed points, at sufficiently low temperatures the phases 
4>s±,4>c- freeze to their minimum positions, and anti- 
symmetric current flows by means of solitons. Again, in 
this case we expect drag resistivity to be proportional to 

p T*/T 



V. ESTIMATE OF PARAMETERS 

The significance of Coulomb drag in a particular exper- 
imental set-up of two ID coupled wires is determined by 
T* and the corresponding minimal length L* = olEq/T* . 
We now estimate these quantities for a double wire sys- 
tem with parameters taken from the experiment by Ya- 
coby et al. Ej], where, using cleaved edge over-growth in 
GaAs-AlGaAs structures, quantum wires of width down 
to linm and length L of order \im were fabricated, with 
adjustable electron density. 

The bare values of the interaction constants gi...,g~i... 
appearing in the fermionic Hamiltonians can be esti- 
mated from the geometry of the experiment. If we as- 
sume the wires to be parallel in a distance d and char- 
acterize the transversal extension of the electron wave- 
function by a length d, intra- and inter-wire potential 
are given by 



V (x) = -(x^ + d 2 )- 1 ' 2 , V(x) = ti -(x 2 + d 2 + d 2 )- 1 ' 2 . 
e e 

The expression for V(x) is approximate on small length- 
scales ~ d. However, since generally kp 1 ~ d, this may 
not lead to large errors in the determination of the 2k p 
scattering parameters, g\,g\. The couplings gi,gi, are 
the Fourier components of V{x) 1 V{x) at q\ = 2kp and 
92,94 -> 0, 



2e< 



2e~ 

K (q z d), gi = K (q t Vd 2 + d 2 ). 
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Taking the limit q\ — > does not raise any problem, 
since the q — ► couplings £/2/4>52/4 appear always as 
differences. Screening by metallic gates used in the ex- 
periment is negligible due to the relative large distance. 
As pointed out by Starykh et al. ( plj , see also Ma- 
han |3^]), the parameter g^ has to be modified for the 
Pauli principle to be satisfied. This modification is how- 
ever small when the 2kp part of the interaction is much 
smaller than the q ~ part. 



Having determined </»..., <fa..., we can calculate the bare 
parameters Ao, K®_, and Ao, K° for the bosonic Hamilto- 
nians by means of the corresponding expressions of Sect. 
Ill and These values are then to be used as initial 
values (x = 0) for the RG equations. 

The system enters the regime of strong coupling when 
A(x) ~ 1 (spin-polarized) or A(ar) ~ 1 (spin-unpolarized). 
Thus, T* and L* are estimated by A(x*) = 1 or X(x*) 
I . respectively. The integration of the RG-equations can 
be done analytically in certain parameter regimes, which 
leads to the expressions (p2|), (|23|), and (|53|) (see appendix 
|A|). Here we also give numerical results. 

One should be aware that this method leads to or- 
der of magnitude estimates rather than precise values 
for T* . Nevertheless, the reliability of this procedure 
can be demonstrated in the special case of spin-polarized 
fermions with K c ^ = 1/2. Then, the bosonic Hamilto- 
nian iJ c _ can be transformed to a Hamiltonian of non- 
interacting, fictitious fermions ( "refermionization" ) [[l6[ , 
which exhibit an energy gap of the order of T* obtained 
by the aforementioned method (see appendix |^) . 

Let us first consider the spin-polarized case. We as- 
sume an inter-wire distance of d — 3d. Table || con- 
tains lg 10 Eq/T* = lg 10 L* jr for wire widths d = ag, 2clb 
and 4as, where as is the effective Bohr radius. With 
as ~ lOnm in GaAs these wire widths are close to the 
experimental values in pl[ | . The TL-parameter of the rel- 
ative charge sector is typically K c - as 0.8 and Ao ~ 0.1. 
The main characteristic is the strong decrease of x* with 
r as long r ~ d, in accordance with the exponential sup- 
pression of T* (see Eq. (|2^)). For larger values of r the 
dependence is rather weak, since here, according to (|23|), 
x* ps ir 3 c(kF)/r s kFd with c(kF)/r s kp w const.. 

The strong coupling regime is not easily accessible. 
Even at very low densities of r _1 = (200n?7i) _1 the 
crossover temperature is still very low: T* ~ 0.01-Eo ~ 
lmK (the corresponding length is L* ~ 10/im). How- 
ever, even when the system is still in the weak coupling 
regime, the drag resistance can be significant. According 
to the Eq. ( |2"c| ) we estimate the drag resistance of a 10^i 
long double wire at T = lOOmAT with d = lOnm, d = 3d, 
and r/a-B = 5,10 and 20 to be of order 0.01, 0.1 and 
lx/i/e 2 , respectively. 

Quantitatively, our findings deviate strongly from the 
estimate for L* given by Nazarov and Averin [ fL3[. For 
parameters d = d = lOnm, D — lOOnm Ref. ]13[] ob- 
tain L* « 0.3/ito, whereas according to our calculations 
for these parameters L* w 1000/zm. The origin of this 
discrepancy of more than three orders of magnitude is 
explained in the next section. 

The results obtained for a double wire of spin- 
unpolarized electrons is shown in Tab. |n] (still d = 3d, 
D = 200nm). For relatively small values r ~ d the 
system scales to strong inter wire couplings A. In this 
regime, x* decreases with increasing r, as in the previous 
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case. A qualitatively different behavior sets in at larger 
r: here, renormalization drives the system towards the 
spin gap fixed point, indicated by a "*" in Tab. |J. Th e 
transition from the former regime to the latter happens 
at r c /a,B = 3.1,4.8, and 7.6 for d/as = 1.0,2.0, and 
4.0, respectively. At these densities, the TL parameter 
assume values A" c _ « 0.5 and K s « 1.3. Since further 
A 0.3, it is unclear, however, to what extent the RG 
equations we use, which are derived for small An, An are 
applicable for this case. 



VI. SUMMARY 

In this paper we used the Tomanaga— Luttinger model 
to analyze Coulomb drag between two ballistic quantum 
wires. We find the drag to be a strong effect, both in its 
magnitude and in its temperature dependence. 

We find that at zero temperature, for all weak 
monotonously decaying repulsive electron-electron inter- 
action the trans-resistance diverges, indicating the forma- 
tion of inter-locked charge density wave ground state in 
the two wires. At low temperature, T <C T* , we predict 
the trans-resistance to depend exponentially on T* /T. 

At high temperature T > T* we predict the trans- 
resistance to show a power law dependence on tempera- 
ture, with the power being determined by the Luttinger 
liquid parameters. For spin polarized electrons, we find 
the power to be 4A C _ — 3, with K c - being the TL- 
parameter corresponding to anti- symmetric charge dis- 
placement. In this, our findings arc in contrast to that 
of Ref. where the exponent was identified with the 
parameter corresponding to symmetric charge displace- 
ment. Quantitatively, the two parameters are different, 
and this difference leads to the big difference in the esti- 
mates for T* . For spin unpolarized electrons, we find the 
power to be 2A C _ — 1. 

The crossover temperature T* depends exponentially 
on parameters. For wires at large separation kpd 3> 1 
it is T* — £ cxp- , where b = 1 for spin- 



An 



2e 2 



-Ko(2fcF d) ~ r s e 



-2k F d 



unpolarized, and 6 = 2 for spin-polarized wires. 
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is exponentially small, as a consequence of which 
K c -(x) « const. = K (Eq. [ll]). With this we obtain 
from d In X/dx = 2 - 2K 

X(x) = X e (2 - 2Ko)x , (Al) 

where x = In E /T. The condition X(x*) w 1 then gives 



T* « E Xo 



1/(2-2X0) 



Eo exp 



kpd 
1-Ko 



(A2) 



For small wire separations (d ~ d) inter and intra 
wire couplings g and g become similar, such that 



U 



(.91 -52+52) 



91 



2irv c - 



2TTVr- 27TU r _ 



An 



This means that the bare couplings Ao, Kq w 1 + Uo lie 
close to the separatrix, but still outside the attractive 
region of the weak coupling fixed points (see Fig. [l]). 
Under renormalization (X(x),K(x)) flows along the sep- 
aratrix towards weaker couplings until the turning point 
at K = 1 is reached. Thenceforward (A, K) flows towards 
the strong coupling regime. 

For an estimate of T* we use the approximate RG 
equations valid for \U\ <C 1, 

dA oxrr dU „ 
— = -2XU, — = -2A 2 . 
dx dx 

These can be easily integrated by use of the constant of 
motion A 2 — U 2 = a 2 , which yields 

X2 — xi = — (arctan Ui/a — arctan JJ-xIa). 
2a 

In particular, we obtain for the "time" x a , at which the 
turning point is reached (x\ — xq = 0,Ui = Uo,x a — 

X2,U 2 =0) 

x a = — arctan Uo/a. 
2a 

In the case under consideration with An rs Uq, the con- 
stant a is a small quantity: 

a 2 «2A (A -C/ )« A 2 ,C/ 2 , 

which will be determined more explicitly later. Because 
of the relative smallness of a we have 



APPENDIX A: 

We begin with the calculation of T* for spin-less elec- 
trons and large wire separation d, such that kpd 3> 1. In 
this case 



" 4a 

For the rise from the turning point to stronger couplings 
A » a, an equally long "time" Xb is needed, so that 

x* w x a + xi, = n/2a 
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and therefore 

T* « Eoe~ n/2a . (A3) 

It remains to calculate a: Taking the expressions 

2e 2 2e 2 
9i = — Ko(gid), ft = — K Q (qid), 

and expanding up to second order in kp d we obtain 



din A 

dx 



= 1 - A' - U s - 2A, 



the last two terms, if Z7° , Ao -C 1 — K®s. Then, for this 
regime, 



din A 

dx 



1 - K 



which gives 



An - C7n = 



1 



-(ffi - 91 +52 - 02 ) 



^f(d/d,\nk F d) (k F d) 2 + O i (k F d), 



X(x) = X e^- K >. 



This leads in the same way as in the corresponding spin- 
less case to a crossover temperature 



whereby 

d - d 2 

f(d/d,\nk F d) = (^lnd/d + 7- l + lnfc F d)(^- - 1). 

d d 2 

(7 ~ 0.577 is Euler's constant) This yields 



E exp - 



2k F d 
1-K9 



a = (2A (A -[/(,)) 



1/2 _ r s i /2 



-^Kl /2 (2k F d)f 1/2 k F d. 



As long A: F J ~ 1, the Bessel-function Ko(2k F d) and / 
are of order one and only logarithmically dependent on 
k F d. Hence in this case 



a = c 1 — s -^k F d 
2tt z 



This result deviates from the former estimation for the 
spin-less case only by a factor 2 in front of k F d. This 
extra factor reflects the fact that here the mean electron 
distance f is exactly half of the one in a spin-less wire, 
when the Fermi momentum is the same. 



APPENDIX B: RELATION TO THE 
MEMORY-FUNCTION FORMALISM 

In this section we elaborate on the relation of the per- 



with c a numerical coefficient of order unity. Inserting turbative calculation of subsection [III D| and a formula 



this into (A3) yields 



E exp - 



r s k F d 



Let us now consider a spin-full double wire with large 
separation d 3> kp 1 . Then again the inter wire backscat- 
tcring coupling is exponentially small, such that also here 



const. 



Due to the smallness of A, the RG equations for the spin- 
sector couplings A, K s decouple from the relative charge 
sector and become identical to the corresponding RG 
equations of a single spin-^ wire: 

A' = -2Z7 S A, 
U' s = -2 A 2 

(for A, \U S \ -C 1). Accordingly, the couplings A, U s with 
bare values A = U° = gi/2irv F scale down to weaker 
coupling: 



for the Coulomb drag resistivity based on the memory- 
function formalism by Zheng and MacDonald (Eq. (12) 
in §§). 

From Eqs. @ and (ph follows that 



d / 8H int {Q) 

9 dn\ 5(f>{x ) 



where M is a constant factor determined by the system 
parameters, and H int is given by Eq. (|l7|). A first order 
expansion in Hi nt — Hi n t(Q) leads to 



/ SH mt \ f°° , /, TT 1 \ 5H int (t = 0), 



1 \ 



where the subscript indicates a thermal average taken 
with respect to Hq. For / oc il — > we can expand Hi nt 
in 

H int {n) « H mt {n = 0) + % [ S -^p\dx, 

V2 J 0(p{x) 



X(x) = U s (x) = 



Ao 



1 + 2X x 



Because of this behavior we can neglect in the differential 
equation for A, 



and obtain 

M f°° 
p = i — = / dt t I dx 

V2Jo 



- SH int (t) SH int (0) . 
5<p(x) ' 8(j>{xo) 



Then, making use of the Kubo-identity 
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Jo dt 



dA 



-ie- f3Ha I d\e XHo ^e- XH ° 
dt 



and partial integration leads to 

V2Jo Jo J \ ty(x ) S(f>(x) / 

Since the functional derivatives of Hi nt are the intra- 
wire force-densities (up to a constant factor), this expres- 
sion essentially equals the corresponding formula (12) of 
Zheng and MacDonald. 



APPENDIX C: REFERMIONIZATION 

Luther and Emery [jl6| observed that in the special case 
K c - = 1/2 the bosonic Hamiltonian i/ c _ can be mapped 
onto an exactly solvable Hamiltonian of fictitious, nonin- 
teracting fermions ("refermionization"). After the canon- 
ical transformation II C _,<^ C _ — ► II' = 2 1 / 2 II C _,<// = 
2~ 1 / 2 (i r _ the Hamiltonian states 



H r _ = — I dx 7T 2 n' 2 

2tt 



' — I dx cos(2<£'), 
2na 



(d x cf>') 

which expressed by fermions 
Or (A;) = 2-1/2 dx Mx)e~ lkx 

i' r = (27ra)~ 1/2 exp(ir(k F x + <p') + iir / dx'W(x')), 

becomes 

H*_ — u c _ rkal(k)a r (k) 



+ -^- V al (k)a-(k - 2k F ) + h, 
2ira z — ' 



This Hamiltonian is diagonziable by a simple Bogoliubov 
transformation. The spectrum 



11, -ki- r l/ (^~) +Uc-(fc-fcF) 2 



exhibits an energy gap of width AE 1 = ,g/7ra at ifcp, 
which coincides up to a factor 2 with the corresponding 
temperature T* as estimated by integrating the RG flow: 
from Eq. (|A§) with A' c _ = 1/2 we get 



2iru c 



9i 
2ira 
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TABLE I. The table lists \g w E /T* (which is also 
lg 10 L*/r) for different values of d/as and r/aB for a 
spin-polarized double wire system 
as = lOnra). 
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FIG. 1 
electrons 
system with 
d< kZ 1 . 



The RG-flow of a double wire system of spin-less 
Point A corresponds to the bare couplings of a 
S> kZ 1 , point B to wires with narrow spacing 




TABLE II. lg 10 Eo/T* = lg 10 L* /a for a spin-unpolarized 
double wire system (d — 3d, D — 200nm, as = lOnm). The 
"*" indicates the zero-drag phase at T = 0. Note that at 
large r/aB the initial values of A, A are not small, and thus 
the validity of the RG equations is questionable. 



FIG. 2. RG flow of the spin-full double wire system in 
the K c — K s = plane (schematic), po and p\ are unstable 
fixed points, p2 is the strong coupling fixed point, and P3 a 
spin-gap fixed point. 
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FIG. 3. Stability of weak coupling fixed points A = A = 0. 
Only fixed points in area / can be stable. In the limit 
A()| <C Ao <C \K S — 1| points in area // flow to the strong 
coupling fixed point P2 (where A — > co) , while points in area 
III or III' flow towards the spin-gap fixed point P3 (where 
A -> 0). 



